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We have investigated the bound states of an electron and positron in superstrong magnetic fields 
typical for neutron stars. The complete relativistic problem of positronium in a strong magnetic 
field has not been succesfully solved up to now. In particular, we have studied the positronium 
when it moves relativistically across the magnetic field. A number of problems which deal with the 
pulsar magnetosphere, as well as the evolution of protoneutron stars, could be considered as a field 
for application. 



I. INTRODUCTION 

Theoretical models for radio-pulsar emission mechanisms include the bound states of a relativistic electron and 
positron in a superstrong magnetic field. The year 1985 was marked by a number of papers, where the authors tried 
to take into account the positronium contribution into the dispersion equation of a photon which propagates in a 
superstrong magnetic field (Leinson & Oraevsky, Q, 0); Herold et al., ||; Shabad & Usov, [Q). Nevertheless, the 
correct solution of the completely relativistic problem for the bound states of an electron and positron in a superstrong 
magnetic field remains unknown. Roller et al. ]5j have formulated the quadratic form of the completely relativistic 
Dirac equation for bound states of an electron and positron in a magnetic field of arbitrary intensity. However, 
they only found a solution for the simplest case, when positronium does not move across the magnetic field. There 
was also the work of Shabad and Usov ||, where they tried to solve the Bethe-Salpeter equation for a positronium 
atom relativistically moving across the superstrong magnetic field. Unfortunately, they did not take into account the 
retardation effect in the interaction between the electron and positron. For this reason, we return to the problem of 
positronium in a superstrong magnetic field. We suggest a solution of the Bethe-Salpeter equation for positronium 
relativistically propagating across the superstrong magnetic field. The results which contradict those given by Shabad 
and Usov are discussed. 

This paper is organized as follows. In Section 2 we state the problem by considering the Bethe-Salpeter equation for 
a bound electron-positron pair interacting with an external magnetic field. The problem can be simplified when the 
adiabatic approximation (discussed in Section 3) holds, since in this case only one Landau level for each particle has 
to be taken into account. In section 4, we concentrate on the particular case when the electron and positron occupy 
the same Landau level. As a particularly important case, we devote Sect. 5 to the ground band of positronium levels, 
and analyze some particular cases. In Appendix A we recall the explicit formulae of the electron (or positron) wave 
function in the presence of an external magnetic field. 



II. BETHE-SALPETER EQUATION FOR A BOUND PAIR OF ELECTRON AND POSITRON IN A 

STRONG MAGNETIC FIELD 

To find the bound states of electron and positron in a strong magnetic field, we start from the Bethe-Salpeter 
equation for the scattering amplitude, to the lowest order in the fine structure's constant e 2 ~ 1/137 

T h (p_, -p+) = ie 2 J y^ Ajtcu (p-,P-) lu T st(jP- - 9, -P+ ~ qht q G q j (-P+, -P+) D MV (q) (1) 

where G (p-,p_^j and G (— P+>— P+) a- re the propagators of the electron and positron, respectively, in an external 
uniform magnetic field, which will be taken as directed along the Z axis. Here and henceforth, a — (+) subscript 



1 



denotes the electron (positron), and summation over repeated indices will be understood. -D M „ (q) , with q = (u>,q) is 
the photon propagator. In Eq. (|l]) p^,p +1 q,p + and p_ are the four-momenta of the interacting particles. 

It is convenient to study this equation in the mixed (coordinates-energy) representation. To this purpose, we define 
the new unknown function 

, , f c£ 3 p_ d 3 p + .. . , . . , . 

Xij (r_, r+ ;£_,-£+) = / — — 3— — 3-cxp (ip_r_) r. y - (p_, -p+)exp ^p+r+) . (2) 

J (27T) (27TJ 

which depend on both individual space coordinates rzp and energies of the electron and the positron. Their energies 
can be written as: 

E E 

£- = —+£, £+ = —-£ (3) 

where the total energy E of the bound pair is an integral of motion in the stationary state we are considering. By 
the use of Fourier transformation of the Green's functions 



G 



(p-,p'_^j = J d 3 r'_d 3 r_G(r_, r'_ ;e_)exp [-ip_r_ + ip'_r'_] (4) 



G (-P+, -P+) = J d 3 r' + d 3 r + G{r' + ,r + ; -e + ) cxp [ip' + r' + - ip+r + ] (5) 

we obtain the following equation: 

(E E \ f 1 1 ' E E 

r_,r+; — + e,-—+ ej = ie 2 J d 3 r'_<2 3 r + G«(r_, r_;— + e)G qj (r' + , r + ; -— + e) 

/d 3 aduo . r , . . N1 „ / , . E E \ 

— — T D„ U ( w,q)exp L*q (r_ - r + )J T^Xrt ^r_,r + ; — + e - w,-— +e - wj 7 t , (6) 

We choose the following gauge : 

A = 0, A = (0, Bx, 0) (7) 
so that B x = By = 0, B z = B. In this way, the Green's function of the electron takes the general form : 

G(r_,r_; £ _) = ^ j ^^^bh (z--Z_) + i k 2 ^-A^lL- (8) 



where E n (k 3 ) = y m| + (fc 3 ) and m„ = m\/l + 2n6, with 6 = B/B . (B = m 2 /e = 4.4 x 10 13 G is the so-called 

critical magnetic field). The electron wave functions "J^nt are bispinors of positive frequency which correspond to 
Landau states of the electron in the magnetic field. Each Landau state of the electron is marked by the number 
n = 0, 1, ... which characterizes the quantized motion across the magnetic field; k^ is the electron momentum along 
the magnetic field; a = ±1 corresponds to the spin projection of the electron S3 = ±1/2. For many applications, 
it is convenient to use wave functions of negative frequency for the states of positron (See Appendix A). Then, the 
positron propagator has the following form: 

r( > r s ^[dk+dkt , , , , ¥-l( X+ - % -fc 3 +)^(.T + - % -fc+) 

G(r,,r, e+)= > / — — exp — ikZiz, — z + )— %kJ [y, — y + )\ — r- — — , (9) 

V+ + +; faj 2tt 2tt ^ l 3V+ +i 2yy+ y+n [ £+ + E n ,{k+) - i0] ' W 



where E n /(k£) = y m^, + (k£) . We assume that the center of mass of the positronium can move relativistically 
in any direction with respect to the magnetic field, i.e., we make no assumption about the center of mass momentum 

P3 = fcf + k+ (10) 

along the magnetic field, but we assume that the longitudinal relative motion of the electron and positron is nonrela- 
tivistic. Therefore, the relative momentum projection 
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k 3 = ^3 — "^L (ii) 

is much smaller than the electron mass. Let us write the electron and positron momenta in the following form : 

fc 3 ~ = ™" P3 + h (12) 

m n + m n i 

k += P 3 -fc 3 (13) 

m n + tyi i 



Thus, we may expand the Landau levels energy as a series of k, 
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" ( 3 ] " ^ x^T^ (Ml, + P|) 3/2 2to " ( } 



where M nn / = m n + m n < . This yields 

2/i 

with the reduced mass 



K(fc3")+^(fc 3 + ) -V / Ml~+^+^ 7 (16) 



_ m n m»> (Ml,+P 3 2 ) 3A - 



III. ADIABATIC APPROXIMATION 



In a sufficiently strong magnetic field i? >> 10 9 G the Larmor radius 



a L = -L= (18) 



is small with respect to the Bohr's radius 

a B = -\ (19) 

because : 

^ = = ^»1. (20) 

When inequality ( pp| ) holds, the Coulomb's binding energy |e| ~ me 4 is much smaller than the distance between 
Landau levels of the electron (or positron). This means that the energies of the bound electron and positron vary 
in a small vicinity of the Landau levels they occupy. In this case, the poles near these Landau levels give the principal 
contribution to the Green's functions of the electron and the positron. Therefore, we can write : 



G n (r_,r ;e_) = V / """ a ~" 2 expftfcT^- - z )+iK(y--y )] — ' ! ^ ■ (21) 

^ [ dk+ dk+ n _ ! ju+rJ ^ iU+fJ „ <7 j(^ - gy-fc 3 ^y(x ± - §, -fc+) 

[e+ + £„,(£;+) -zO] 

(22) 



G™'(r + ,r + ;e + ) =J] / ^ cxpHfcj-(z+ - z + ) - ik+{y' + - y+)] 
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Generally speaking, in the case of n ^ n' two terms corresponding to the same total energy E — e + + e_ give the 
major contribution to the product of Green functions in the right-hand side of the Bethe-Salpeter equation. The first 
of them corresponds to an electron occupying the Landau level n, while the positron occupies the level n . The second 
term is the state with the same energy, but with the electron occupying the level n and the positron in the state n. 
Within this adiabatic approximation, it is possible to characterize the positronium states in a strong magnetic field 
by the almost-good quantum numbers n and n , and by other additional constants, which characterize the motion of 
the center of mass and the relative motion of the electron and positron along the magnetic field . 



IV. EQUATION FOR THE BOUND STATE WAVE FUNCTION FOR EQUAL LANDAU NUMBERS 

If both the electron and the positron occupy the same Landau level, then virtual photons with small values of 
~ me <C 
propagator 



ui ~ me 4 -C |q| ~ me 2 give the major contribution to the e + e interaction. Therefore, we can neglect u> in the photon 



47T 4-7T 

Dfj, u ( u, q) = g^—^ J - ( 23 ) 

lo z — <\ Z q 

This yields, for n = n'\ 

Xij (r-,r+; | + e, ~ + e^j = -ie 2 J d 3 r'_d 3 r' + G^(r_,r_;j + e)G n qj (v' + ,v + - ~ + e) 

X / q^ 6XP ^ ~ r + )] l S Xst ( r -' I + W ' ~f + 7t < (24) 

By performing integration over de/2n, we obtain the wave function of the bound pair 

f de ( E E \ 

$(r_,r+;£) = / — Xij r_,r+; - + e, -- + £ (25) 



which verifies the following equation: 



/tflr rfir' dp W, R 

| r> -_ r/+ + gG"(r-,r_;-+ £ ) 7 ^(rL,rV;^) 7 ,G"(rV,r + ;--+e) (26) 

Integration over de /2tt in the right-hand side can be done with the help of the following formula : 

d£ 1 1 (27) 



Loo 2tt (e + E/2 - E n (kz) + iO) (e - E/2 + E n ,{k+) - ffl) E- E n (k^) - E n ,(k 
Thus, one has 

QTq 3 = J ^G3(r_,r'_;|+ £ ')G«;(r + ,r; ; -| + e ')=^^i J 



ds „„ , i E i, i , i E i , x -> . /" dk 3 dk 2 dk^ dk^ 

~2tt 2~k 2~7r 2~k~ 

exp[ik~(z_ - z_)+ik^(y- - y'_) - ifcj~(>+ - z+) - ife^(y+ - y+)] 



E-E n (hz)-E n ,(k+) 

x^+\x. + ^)^\ x + ^)^r)(x'-^w { r)(x + -% (28) 

eS y nff 1 " eS y + " CT V + eB J K ' 

Since the relative motion of the bound pair along the magnetic field is nonrelativistic, we neglect ~ me 2 inside 
the bispinors. The Landau level energies of the electron and positron are given by Eqs. ( |l4] , |l5| ). 
Let us introduce the new coordinates : 

„ m„z— + m„i z+ 1, . 1, 

Z= n , " + , K =-(»-+!/+), X= 
m„ + m n » 2 z 

z = y = y_-y + , x = x_-a; + (29) 
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and new momenta 



m„ z + m„ i z i l,i i . 1 . / / , 

Zi = „ + ^ = y _ + ^ = , + 

m„ + m„' 2 2 

z\ = z_ - z + , yi=y_-y + , x\=x_-x + (30) 



P 3 = k 3 + k+, P 2 = k 2 + k+ 

h = m " fc3 '' m " fc \ k 2 = i(fc 2 - - kt) (31) 
m„ + m< 2 



The expression for Q nn i reads now : 



. dP 3 dk 3 dP 2 dk 2 exp{iP 3 (Z - gi) + iP 2 (F - Yi) + zA: 3 (z ~ zi) + ik 2 {y - yi)} 
^ u <" 1 2tt 2tt 2tt 2tt c _ _ 



2/i„„/ 2M „„' 

a; , P 2 + 2/c 2 ,-ri(+) , v , a?i , P 2 + 2k 2 



x**-?^ - 5i - P2 - 2fc2 )¥ J( ,i(X - £ - P2 ~ 2k2 ) (32) 

" ff v 2 2eB ' n a y 2 2eB J y ' 

and the positronium wave function obeys the following equation: 

At t> ip\ 2 f dZ 1 dY 1 dX 1 d 3 r 1 dP 3 dk 3 dP 2 dk 2 e X p{iP 3 (Z-Z 1 )+iP 2 (Y-Y 1 )} . 
Z(r,K;E)=e J _____ g g exp{*fc 3 (* " *i) +Mv ~ »)} 

x y y {x + x - + P2+2k2 w l ^ ] (x 1 + ^ + P2+2fc2 ) 

_ ^_ n CT v ^ 2 2eP 2 2eP- ; 

" a' 

xy$(r 1> R 1 ;_)7,*£7 ) (^ - f - ^J^)^(^ " f " ^J^)- (33) 

To separate the motion of the center of mass we subtract a bilocal phase (Avron et al., [fij) 

PR+erA(R) (34) 

from <i>, with the help of the gauge choice Eq. (Q). Thus, one has 

$(r_,r+;P) = cxp[iP$Z + iP$Y - ieBXy]<f>(r) (35) 

By using this expression one can perform the integration on the right-hand side of Eq. ([33]) over dZ\dY\. Thus, we 
obtain the S- functions S(P 3 — P®)5(P 2 — P§ ), which allow for integration over dP 3 dP 2 . Therefore, we will substitute 
everywhere the constants of motion P® , P 2 instead of P 3 and P 2 . We consider the reference frame P3 = 0. Thus, the 
equation for the stationary states becomes the following: 

Ti P - / -^—exp(ik 2 y + ieBXy) / dX x exp (~ik 2 y 1 - ieBX x y\) exp{ik 3 (z - Zi)} 

r\ 2n £ b_ J 2n J 

x W, T ,(+) r v- | * | P " + 2fc 2 X! P 2 ° + 2fc 2 



f-1, a;i P? - 2fc 2 ,—(-), x P 9 °~2fc 2x 



The right-hand side of Eq.(p6[) depends on z only through the exponential exp(ifc3z). The action of the operator 

1 8 2 

— -r^- £ - (37) 

on Eq.(|3^) cancels the denominator on the right-hand side. Subsequent integration over dk 3 yields S(z — z\), and 
integration over dz\ becomes trivial. Finally, we obtain : 
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/dx dy f dk f 

12 1 '2 =f / "2^ exp ^ 2y + ieBX y) / ^1 CX P (- ik 2yi - ieBX^i) 
\j x-^ -\- y-y z J J 



2 ^ 2 )*^(*i + ^ + P2+2fc2 ' 



2eB 



x 7 ^(r 1 ) 7 ^rj(^i- y 
Let us now introduce the new integration variables 



(-) (x *i P 2-^ 2^{-) tv x P°-2k 2 , 

Vna'V 1 2 



2eB 



t = X + 



h2_ 
eB' 



\ = X l 



eB 



(38) 



(39) 



and denote by xq the distance between the orbiting centers of the electron and positron in the plane orthogonal to 
the magnetic field: 



X 



Pi 

eB 



(40) 



This distance depends on the transverse momentum P® > which characterizes the motion of the mass center across the 
magnetic field. Then 



I dx 1 dy 1 
J y/xj +yl + z 



dAex P (-i e BAy 1 )*it ) (A+^^)7^(a ;i , 2/1 ,^<-J(A-^±^) 



As follows from the last equation, the wave function describing the relative motion has the form : 



^ - E E (*) / § cx p ( jeS ^) k 



x + Zq — (-) _ x + zp , 

o / ncr' VS r> 



(41) 



(42) 



where the functions /™ CT , depend only on the relative coordinate z. These functions are solutions of the following set 
of Schrodingcr-likc equations: 



1 d 2 



(z, x ) + E V$ (*. x °)f™> (z, x ) = eft,,, (z, x ) 



(43) 



Here, the effective potentials represent the Coulomb interaction, averaged over the fast transverse motion of the 
interacting particles : 

/ , n _ 2 eB [ dxdyd£,d\ exp [ieB (£ - A) y] 

V aa' \ Z ) — e 



2ir 



\J {x — Xq) 2 +y 2 + Z 2 



xTr 



(44) 



In general, the matrix V££i (z) consists on 16 elements corresponding to different spin orientations of the two particles. 
However, in the adiabatic approximation, most of them are zero. After integration over d\d£, we obtain the effective 
potentials in the following form: 



dxdy exp(— p 2 /2) 

y/(x- X ) 2 + y 2 + Z 2 



(45) 



where 



Rli (x, y) = cxp (p 2 /2) J ^ exp [teB(t X)y] 



(46) 



with p = ^/e£> (x 2 + y 2 ). 
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V. THE GROUND BAND OF POSITRONIUM STATES IN A SUPERSTRONG MAGNETIC FIELD. 



The lowest band of positronium levels corresponds to the case when both the electron and the positron occupy 
the ground Landau level with n = n' = 0. Only a spin combination, spin-down for the electron and spin-up for the 
positron, is possible in this case. A simple calculation, using the wave functions of Appendix A, gives the following 
result: 

R}}(x,y) = R\\ (x, y)=R\\(x, y) = R\\ [x, y) = -R\\(x, y) = -R\\ (x, y) = R}\ (x, y) = (47) 

and 

R[](x,y) = l. (48) 

Consequently, the wave functions of the lowest band, which will be labeled by the quantum number I, have the 
following form : 

*g(r) = /<?•«>(*, * ) / le^CieB^Wft+i^J^tt-^) (49) 

Using the formulae for bispinors corresponding to Landau states (see Appendix A), we can integrate this expression 
over d£. The result is : 

<&) = f^\z,x ) ] ^exp{-^[(x + x ) 2 + y 2 ]} J( 1 ) (50) 

The wave function for the relative motion along the magnetic field can be found by solving the following Schrodinger 
equation : 

- 2/ ! 00 dz 2 ffi ^' + V $ ^' X °W§' 1) ^' x °* =e oo( x o)fi°' l) ( z ' x o) ( 51 ) 
with fiQQ — to/2. The effective potential is 




_2e2 />« pdpcxp(-p 2 /2) I . / 4px a L 
li 1 -""- n J Q y/(pa L -x Q )* + z* TV (pa L -x o y+zi 



V 00 (z, xo) = V$(z, mi — / ' ' , -'„ K j , , ^± r | ( r,i> i 

and 



/■7r/2 in 

K(fc)=/ - ^ 2 (53) 
Jo J I - k 2 sin 2 9 



■"i 2 de 

K (fc) = ; 
Jo 

is the elliptic integral of the first kind. 

Eq. ( p2| ) depends on the parameter xo = P® /eB, therefore the eigen-functions /j? , as well as the discrete spectrum 

eigen-energies e l 00 , depend on the quantity P^ , which characterizes the relativistic motion of the center of mass across 
the magnetic field. Unfortunately the integration in Eq. ( p2[ ) can not be done analytically. By this reason, we consider 
several limiting cases. 



A. The case of small xo . 



We first consider the case of a small distance between the centers of Landau orbits in the plane orthogonal to the 
magnetic field: 



' eBxn 



s/eB 



< 1 



(54) 



This yields: 
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Voo(z) 



neB 




exp 



(55) 



with 



erf (x) 



eft exp (—t 



(56) 



being the error function. At distances \z\ 3> ol, the potential behaves as the one-dimensional Coulomb potential 



while for \z\ <C eti, goes to a constant 



Vqq(z) ps - — 



(57) 



K,o(0) = -Jf- 



(58) 



(This result differs from that given by Usov and Shabad by the factor . 

For shallow levels, the characteristic size of the atom along the Oz axis is of the order as 3> cll , so that the 
potential can be considered to have the form Eq. (p7|), and the energy levels should be close to Balmer's spectrum. 
This statement is not valid for deep levels. Eq. (pjl), with the potential Eq. (|55|), can be solved numerically. However, 
an analytical estimate is of interest. To obtain such an estimate, we replace the potential Eq. ( [sBl) by a simple 
function of the form 



V Q0 (z) 



(59) 



which slightly differs from Eq. ( |55|) only in the region \z\ ~ a^, and has the same asymptotic forms. 



1. Discrete spectrum. 

Solutions to the Schrodinger equation with this potential have been investigated by Loudon Q for negative energy. 
For completeness, we quote shortly this calculations, which will be used later for investigation of the states of positive 
energy. We define, as in H, a dimensionless quantity a by writing 



1 



with aoo = (/^ooe 2 ) 1 = 2<Zb, and replace the independent variable z in Eq. ( |5l| ) by 

2 ( [2 \ 

cll + z , for z > 



(60) 



aa 00 \ V 7r 

z' = — (\/-a L -z|, for z < (61) 

aa 00 y V 7T / 

whereupon, for z 7^ 0, the equation takes the Whittaker's form of the confluent hypergeometric equation 

This equation has two independent solutions, given by Whittaker's functions. The first of them 

M a 1 (z') = exp (-z'/2) z'$(l - a, 2; z) (63) 
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where $(1 — a, 2; z') is the confluent hypergeometric function of the first kind, diverges as (z')~ a exp (z'/2) for large 
z. Since for a bound state any solution of this equation must go to zero when \z'\ tends to infinity, we choose the 
second solution, which for z > is 



W a> i (V) = exp (-z'/2) z'tf (1 - a, 2; z') 
where ^(1 — a, 2; x) is the confluent hypergeometric function of the second kind : 

*(1 - a, 2; x) = — — { $(1 - a, 2;x)\l gx + ^ (1 - a) - $ (1) - V (2)] - — + V 1 " ' 

I ( — a) rim 



ax ^-^ r\ (r + 1)! 



(64) 



A/} (65) 



Here 



n + 1 — a n+1 n + 2 



(c) r = 



r(c + r) 

r( c ) 



(66) 



The function ip (s) = V (s) /T (s) is the logarithmic derivative of the gamma function. The solutions for positive and 
negative z can be joined together to form either even or odd wave-functions. For an odd state we require 



2 2a L 
7r aa 00 



(67) 



while for an even state 



dz 



(68) 



To find solutions to Eqs. (|67| ) and (|6§|), which give the eigenvalues of the system, we keep only terms which are 
dominant when 



2 2a L 
7r aa o 

is very small, and a is close to a positive integer. The eigenvalue conditions then become : 
Odd state: 

*(l-a)JI**.-i = 



Even state: 



aa 00 V 7T 

Assuming a — > I where I = 1, 2, ... we can replace 



7r aaoo a 



2 (2 

log T — \j—a L + tp(l-a) = 



(69) 



(70) 



(71) 



ip(l-a) 



Then, the quantum defects Si — a — I are given by 
Odd state: 



Even state: 



a — I 



jr 2 2 
a 00 V 7T 



(72) 



(73) 



i 2 2 
log t \ -a L 

la 00 V 7T 



(74) 
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For the adiabatic case one has Si <C 1 if log(a o/aL) ^> 1. 

In addition to this series of states having their quantum numbers a close to positive integers, there is another state 
having a close to zero. For such a value of a, ip (1 — a) is no longer an important term in Eq. (fn|), and 1/q becomes 
the dominant term in a. In this case, the equation for the eigenvalues becomes 



log- 



aa 00 V 7T 



2 1 
-a L = - — 



2a 



(75) 



The quantum defect Si = a for the state Z = can be found from Eq. ([75]) by iteration. To the first order in a^/aoo 
its binding energy is given by 



jo) _ JO) 



me log t / 

1 W 2 2a,, 



7T a o 



(76) 



To this state there corresponds the even wave- function 



In 



a oo 



7T O00 

2~o~l~ 



-.1/2 



exp 



In 



7T aoo 
2"a7 



(77) 



2. Continuous spectrum. 

The spectrum for positive energies is continuous. Now, the variable z', as well as the a value defined by Eq. (|6C 
are imaginary quantities. Let us write them as 



mi 



(78) 



with k = ^/2/ioo£j and define: 



z' = —2ik \ —cll + z 



for z > 



z 1 = 2ik \ —aL 



for z < 



(79) 



In the region z' > 0, Eq. (|62|) has two linearly-independent solutions, given by the Whittaker's functions 



f (fc), 





-2ifc 















(80) 



/f (^) = C 2 M « , 



-2ik 



2 

—a L + z 

IT 



(81) 



where Ci and C2 are constants, which should be determined by the normalization conditions. To find these constants, 
let us consider the asymptotic behavior of Whittaker's functions when z — > 00 



W i 1 

fca 00 ' 2 



-2ik 



-cll + z 
ir 



exp 



3tt 
2fca o 



. i 2 

exp ikz H In 2kz + i\ —kciL 

\ ka 00 V 7T 



(82) 



M < 1 -2ik 



-a L + z 

7T 



y fca 00 J 



exp 



-ifcz 



ka Q0 



\n2kz 



2 , 3tt 
-ka L + 1 — 

7T Z 



iSk 



(83) 
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where 



ka m 



5 k = arg T 

If we choose to normalize the eigen-functions in the following way : 

/DO 
f( k \z)f (k '\z)dz = 2ir5{k-k') 
-oo 

then the normalization factors are 



Cl ( 2fca 00 ) 



C 2 



1 



fca 00 



fcaoo 



exp 



3tt \ 
2ka Q0 ) 



(84) 



(85) 



(86) 



Indeed, the asymptotic form of the wave functions in this case are in a good agreement with the general form of 
the normalized one-dimensional wave functions for a continuous spectrum, in the form exp(±ifcz). The logarithmic 
term in the exponentially grows much slower than z. This term is not important when calculating the normalization 
integral, which diverges at infinity. 



B. The case of large xo. 

The bound pair of an electron and positron with a small magnitude of P2 in a superstrong magnetic field is only of 
academic interest. The known mechanism of positronium production in a pulsar magnetosphere @, [@], f§ assumes 
that the created bound pair has a quasi momentum P® > 2m. For a typical pulsar magnetic field B < 0.1i?o, this 
corresponds to the opposite limit : 

pQ 

eBx = -j^= > 1 (87) 
veB 



In this limit, the effective potential Eq. ([32|) takes the simple analytical form 

V 00 (z/xo) « -—-=L= (88) 
xo ^l + z 2 /xl 

For P® > 2m and B = O.IBq, i.e. \J eBxo > 6, the curve given by Eq.(jS8|) practically coincides with the exact 
potential. In Fig. 1 we show the potential, in units of me 4 , versus the dimensionless distance me 2 z, for two values of 
P® ■ We also show, for comparison, the approximation 

\z\ + x 

used by Usov and Shabad. 

In spite of the relatively simple expression given by Eq. (|S8|), there are not available analytical solutions of the 
Schrodinger problem for stationary states in this potential. For this reason, we have made a numerical calculation in 
order to determine the energy e l Q0 for bound states. The results, for the first five levels, are shown in Table 1, where 
three values of the dimensionless parameter £ , defined as 

£o = (90) 

were considered. For large values £0 3> 1, these results can be approximated by the following formula: 

. 1/ - l.r - s;„ - !/ -!- 1)'" 

'00 



V(Z + l) 2 + 8£ -(Z + l) 

mi 



(91) 
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VI. CONCLUSIONS 



We have investigated the bound states of an electron and positron in the presence of a superstrong magnetic field 
(B >> 10 9 G), as it can appear in some neutron stars. We have given a completely relativistic description of the 
positronium motion across the magnetic field. Our starting point is the Bethe-Salpeter equation for the positron and 
the electron at the lowest order in the electromagnetic coupling constant. The effects of the strong external magnetic 
field are incorporated through the exact solutions of the Dirac equation for the interacting particles. The Bethe- 
Salpeter equation then involves a summation over all possible Landau levels of those particles. As we have shown, 
however, this equation can be transformed into a set of coupled Schrodinger-like equations under the hypothesis of 
the so-called adiabatic approximation, which is valid for superstrong magnetic fields. In this case, only one Landau 
level per interacting particle is relevant. 

We have concentrated ourselves to some particular cases of particular interest, like the ground band of positronium, 
where we found some differences with previous results. 
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APPENDIX A: SOLUTIONS OF THE DIRAC EQUATION IN A CONSTANT MAGNETIC FIELD 

The wave functions we use are solutions of the Dirac equation in a constant magnetic field directed along the z-axis. 
By the use of the asymmetric (Landau) gauge Eq. ([?]), these four wave- functions, with definite third-component spin 
direction, can be expressed in terms of the stationary wave functions as: 



v 7 L2L3 



exp(-ieE n t) exp [it (k 3 z + k 2 y)} ^(Ce, k 3 ) 



(Al) 



The wave functions are normalized in a volume (L1L2L3) , n is the quantum number of the Landau level, a =],[ 
is the spin projection along the z-axis, e = +1 (e = —1) indicates the electron (positron) states, and k 2 ,k 3 are the 
momenta in the y and z directions, respectively. The functions ^/^(^e, fc 3 ), with 



'eB x ± 



•D 



(A2) 



are given by 



"I 



/-ik 3y /E° - my n _i( ^ + ) 
\-k 3y /E°+m<p n (£ + ) 



(A3) 



>*1 



Cl' 2 



( i E n + £°)y/fl° +W»-l(e+) \ 

-ik 3y /E° - m(p n {£_+) 

k 3y /E°+m< p n - 1 (S + ) 
\i(E n + E°)^Ej^; i p n (Z + ) ) 



(A4) 



/faV^S + myn-xte-) 



C l/2 



-i(E n +E° )yE°- mtp n (S-) 
(E n + E°)y/E°+rrup n -i(Z-) 



J 



(A5) 



nl 



c l/2 



-h^Ej 



mip n - 



\(E n +E° n )^EfT^MC- 



) / 



(A6) 
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with 



Cn =4E n E°(E n+ E°) (A?) 



E n = ^m 2 + 2neB + k 2 (A8) 



E° n = ^Jm 2 + 2neB = m n (A9) 
The functions <p n {Q are the eigenfunctions of the one-dimensional harmonic oscillator, normalized with respect to x. 
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TABLE I. Bounding energy of the first positronium energy levels, for different values of the parameter £o 





I = 


I = 1 


I = 2 


I = 3 


I = 4 


0.1 


-1.9093 


-0.2371 


-0.1283 


-0.0609 


-0.0433 


1 


-0.3349 


-0.1374 


-0.0757 


-0.0463 


-0.0318 


10 


-0.0432 


-0.0319 


-0.0240 


-0.0186 


-0.0147 




-2-10 12 



2 

me z 

FIG. 1. Interaction potential of positronium, in units of me 4 , as given by our Eq. ( |S8| ) (solid lines). We also show for 
comparison the approximation, Eq. ( p9| ) used by Usov and Shabad (dotted lines). In both cases, the deepest curve corresponds 
to = 2m, and the other one to i2 = 4m. 
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